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In this contribution an improved analysis is described to extract precise charm and 
bottom quark masses from experimental and theoretical moments of the photon polar- 
ization function. The obtained MS mass values read m c (3 GeV) = 0.986(13) GeV and 
m b (W GeV) = 3.609(25) GeV. 

1 Introduction 

The theory of strong interaction has the strong coupling constant and the quark masses as 
fundamental input parameters. The latter constitute an essential input for the evaluation 
of weak decay rates of heavy mesons and for quarkonium spectroscopy. Furthermore, decay 
rates and branching ratios of a light Higgs boson — as suggested by electroweak precision 
measurements — depend critically on the masses of the charm and bottom quarks. Last 
not least, confronting the predictions for these masses with experiment is an important task 
for all variants of Grand Unified Theories. To deduce the values in a consistent way from 
different experimental investigations and with utmost precision is thus a must for current 
phenomenology. 

The method described in this contribution goes back to 1977 Q] and was applied to 
next-to-next-to-leading order (NNLO) in Ref. [2]. The NNNLO analysis, including updated 
experimental input, was presented in Ref. [3]. 

2 Moments 

The basic object which enters our analysis is the photon polarization function defined 
through 

+ n(9 2 ) = i J dxe^(0\T JtI (x)ji(0)\0), (1) 

with being the electromagnetic current. The normalized total cross section for hadron 
production in e + e" annihilation is then given by 

„, „ a(e + e~ — > hadrons) r „, Nl 

R(s) = -i - = 127rlm[n(g 2 = s + ie)] , (2) 

Opt 

where a pt = 47ra 2 /(3s). In the following we add a subscript Q to indicate the contribution 
from the heavy quark Q. 

The idea for extracting a quark mass value rriQ is based on moments constructed from 
ITq. On one hand one can compute the Taylor expansion of ng(q 2 ) around q 2 = and 
obtain the so-called "theory-moments" from 
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The three- loop contribution to TlQ(q 2 ) up to n = 8 within QCD has been computed in 
Refs. [H [5] and the four-loop calculation for n — and n = 1 has been performed in 
Refs. [H[7]. In the analysis of Ref. [3] also two-loop QED corrections and non-perturbative 
contributions have been considered. The latter shows a visible effect only in the case of the 
charm quark. 

i 

From dimensional considerations wc have tuq ~ (A4 n ) 2n which implies a stronger de- 
pendence of rriQ on variations of M. n for smaller values of n. Furthermore, higher values 
of n require a careful theoretical treatment of the threshold region and the construction of 
an effective theory. The analysis performed in Ref. [3] is restricted to n = 1,2, 3 and 4. 
Note that precise mass values can only be obtained for the three lowest moments since the 
non-perturbative contributions become too big already for n = 4. 

One of the major advantages of the method discussed in this paper is that we can 
adopt the MS scheme for the quark mass entering Eq. ([3]) and thus directly extract the 
corresponding value for the mass. 

In order to extract experimental moments one exploits the analyticity of Hq and arrives 

at 

Mn = J^RqW, (4) 

where Rq naturally divides into three parts: At lower energies one has the narrow resonances 
which are the J/'F and for charm the T(nS) (n = 1, ... ,4) in the case of the bottom 
quark. The corresponding contributions to M. n are obtained with the help of the narrow 
width approximation for R(s) 

R (s) = 6(s -M R ), (5) 

where the electronic widths T ee are known at the 1-2% level. 

The second part is called threshold region and extends in the case of the charm quark 
from 3.73 GeV to about 5 GeV. In this region the cross section shows a rapid variation and 
can not be described by perturbation theory. Measurements from the BES collaboration 
from 2001 [5] and 2006 [9] provide excellent data for R(s) with an uncertainty of about 4%. 
In order to obtain R c one has to subtract the contribution from the light quarks which is 
explained in detail in Ref. [3]. 

The treatment of the bottom threshold region is quite similar. Measurements of R from 
threshold up to 11.24 GeV have been performed by the CLEO Collaboration more than 
20 years ago [TU], with a systematic error of 6%. No radiative corrections were applied. 
The average value derived from the four data points below threshold amounts to R = 
4.559±0.034(stat.) which is 28% larger than the prediction from perturbative QCD (pQCD). 
However, a later result of CLEO [IT] at practically the same energy, i?(10.52 GeV) = 
3.56±0.01±0.07, is significantly more precise and in perfect agreement with theory. Applying 
a rescaling factor of 1/1.28 to the old CLEO data not only enforces agreement between old 
and new CLEO data and pQCD in the region below the T(45), it leads, in addition, also to 
excellent agreement between theory and experiment above threshold around 11.2 GeV where 
pQCD should be applicable also to bottom production. Further support to our approach 
is provided by the CLEO measurement of the cross section for bottom quark production at 

= 10.865 GeV which is given by a b {*/s = 10.865 GeV) = 0.301 ± 0.002 ± 0.039 nb [12]. 
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Figure 1: In the left plot the data for R(s) are shown as published in Refs. [10] (circles) 
and [TTJ (triangle). The black curves are the predictions from pQCD outside the resonance 
region. In the right plot the older data from [TU] are rescaled by a factor 1/1.28. 



The central value can be converted to i?t (10.865 GeV) = 0.409. On the other hand, if 
one extracts i?h(10.865 GeV) from the rescaled CLEO data from 1984 [TUj one obtains 
Rb (10.865 GeV) = 0.425 which deviates by less than 4% from the recent result [12] . 

In Fig. [JJ the original and the rescaled data from [TU] is shown and compared to pQCD 
and data point from [11] . We thus extract the threshold contribution to the moments from 
the interval 10.62 GeV < y/s < 11.24 GeV by applying the rescaling factor to the data, 
subtract the "background" from u, d, s and c quarks and attribute a systematic error of 
10% to the result. 

The third contribution to the experimental moment is provided by the so-called contin- 
uum region which for the charm and bottom quark starts above 4.8 GeV and 11.24 GeV, 
respectively. In both cases there is no precise experimental data available. On the other 
hand, pQCD is supposed to work very well in these energy regions, in particular since Rq(s) 
is known to order a 2 s including the full quark mass dependence and to order including 
quartic mass effects. For recent compilations we refer to Refs. [T3] [j"4l [15] and would like 
to mention the Fortran program rhad [15] which provides a convenient platform to access 
easily the various radiative corrections. 



3 Quark masses 

Equating the theoretical and experimental moments of Eqs. ([3]) and (f4]), adopting yu = 3 GeV 
(/i = 10 GeV) for the charm (bottom) quark and solving for the quark mass leads to the 
results which are shown in Fig. [2] in graphical form0 It is nicely seen that the results for 
ttlq further stabilize when going from three to four loops. At the same time the uncertainty 
is considerably reduced. Furthermore, the preference for the first three moments is clearly 
visible. Also the analysis for n — 2 and n — 3 leads to small errors, even if we include the 
uncertainty from the yet uncalculated four-loop contributions]^) We emphasize the remark- 
able consistency between the three results which we consider as additional confirmation of 
our approach. 

a The numerical results including a detailed error analysis can be found in Ref. \3\. 
b See Ref. [3] for details. 
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Figure 2: m c (3 GcV) (left) and mi, (10 GeV) (right) for n = 1,2,3 and 4. For each value 
of n the results from left to right correspond the inclusion of terms of order oP a , a\, a 2 s and 
oPg to the theory-moments. Note, that for n = 3 and n = 4 the uncertainties can not be 
determined in those cases where only the two-loop corrections of order a s are included into 
the coefficients C n as the equation cannot be solved for the quark mass. 



The final result for the MS-masses read m c (3 GeV) = 0.986(13) GeV and m&(10 GeV) = 
3.609(25) GeV. They can be translated into m c (m c ) = 1.286(13) GeV and mb(mb) = 
4.164(25) GeV. This analysis is consistent with but significantly more precise than a similar 
previous study. 
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